Interaction-induced localization of fermionic mobile impurities in a 

Larkin-Ovchinnikov superfluid 
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We theoretically investigate the interplay between the fermionic mobile impurity atoms and a 
Larkin-Ovchinnikov (LO) superfluid in a two dimensional optical lattice. We find that the impurity 
atoms get localized and can form pairs when the interaction between the impurity atoms and 
the LO superfluid is strong enough. These features are due to the phenomena of self-localization 
whose underlying mechanism is revealed by an effective model. The impurity atoms with finite 
concentrations can drive the transition from a two-dimensional-checkerboard-like LO state to a 
quasi-one-dimensional-stripe-like one. Experimental preparations to observe these features are also 
discussed. 
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Introduction-lurpniiiies or defects which are unavoid- 
able in the real materials always affect the physical prop- 
erties of these systems in different kinds of ways. The 
impurity problems are particularly important in the con- 
text of superconductors where investigations on the im- 
purity effects can help us to understand the experimental 
anomalies, probe the pairing symmetries and reveal the 
competing electronic correlations [1]. The impurities or 
defects can also be introduced into the cold atomic sys- 
tems either by establishing a fine-grained optical speckle 
field[2] or trapping impurity atoms[3-5]. Compared to 
the condensed matter system, these impurities or de- 
fects are easier to be controlled in a cold atomic sys- 
tem. The experimental achievement has attracted in- 
creasing theoretical interest on the impurity effects in 
a two-component superfluid Fermi gas [6]. Despite the 
impurity effects in a superfluid Fermi gas share some 
common features with that in the superconductors, they 
could have their unique properties. Especially, consider- 
ing the impurity atoms trapped in a Fermi gas are mobile, 
the effect induced by them is actually an interactive pro- 
cess which involves the interactions between the impurity 
atoms and the superfluid Fermi gas, making impressions 
on both of them. 

Such an interactive process is particular interesting 
when the two component fermions have different popula- 
tions. Specifically, by tuning the population imbalance in 
a superfluid Fermi gas 0-01 1 it has been found in exper- 
iments that the magnetization or population imbalance 
can coexist with the superfluid, indicating the possible 
existence of the Fulde-Ferrell-Larkin- O vchinnikov [lo[ 11 1 



state. Since it has been well established that the FFLO 
state is more sensitive to nonmagnetic impurities than 
conventional s-wave superconductors or superfluid fl^, 
the mobile impurity atoms are more likely to induce ex- 
otic physical consequences in the FFLO state. 

Here we focus on the interplay between the fermionic 
mobile impurity atoms and a Larkin-Ovchinnikov(LO) 



superfluid in a two-dimensional optical lattice. We find 
that there exits a critical value of the interaction be- 
tween the impurity and the LO superfiuid, above which 
the extended impurity states become localized. The un- 
derlying mechanism of the localization can be captured 
by our proposed effective model. Paired states of impu- 
rity atoms can be formed by occupying the bonding and 
anti-bonding states simultaneously when the number of 
impurities is slightly larger than the saddle points of the 
LO superfluid. On the other hand, novel structures of 
the LO superfluid can be induced and a transition from 
a two-dimensional-checkerboard pattern of the order pa- 
rameters to a quasi-one-dimensional-stripe one can be 
driven by flnite concentration of impurity atoms. This 
provides an indirect method to modify the patten of the 
LO superfluid experimentally. 

Model Hamiltonian-We consider a two-component su- 
perfluid Fermi gas with mobile impurities in a two di- 
mensional optical lattice, which can be described by the 
following Hamiltonian: 



Hcd 



(1) 



H. 



td,i]d\dj and Hcd = f^cd I]^ • He 



d — ^ X/<ij"> 

describes a superfluid Fermi gas, in which ^ represents 
the creation operator of the c-atom with spin-cr at site i 
and Uc is the attractive interaction between c-atoms. tc 
is the hopping integral of c-atoms between two nearest 
neighboring sites. Hd describes the femionic mobile im- 
purity atoms in which td is the hopping integral between 
the nearest neighboring sites. H^d is the interaction be- 
tween the c-atoms and the impurities. Here we consider 
the repulsive interaction with Ucd > 0. 

We decouple model (1) by introducing mean-fleld order 
parameters = -Uc{ci^iCi^^), nl^ = (4 a'^^'y) a-nd nf = 
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Figure 1. (Color online) Contour-plot for (a): A,; without the 
impurities, (b) and (c): density distribution of five impurity 
atoms with Ucd ~ 1-0 and Ucd = 4.0, respectively. Here 
U = 0.3 and f/c = 4.0. 



Then we have equation for c atoms: 



H: 



A* 



(3) 



and impurity atoms: '^jHf^Wj^n — E^Wi^n, where 



H: 



H: 



■d ^ 
ij 
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„ . Ay = AiSij and 

-td.ij + Ucdn'idij. The particle number per site and 



the on-site s-wave pairing order parameter are given by: 



hi 



E 

A» = -UcJ2n^^,nV*„&i-E^) and nf = J2n=i 
where Q{x) is the step function and the magnetization 

Nd = J2i^i is the total 



is defined as m, 



number of the fermionic impurity atoms which occupy 
the lowest Nd energy levels. In our calculation, we choose 
a lattice size of 32 x 32 with periodic boundary condition 
and use = 1 as the energy unit and lattice constant 
a = 1 as the length unit. The total particle number 
and population imbalance of the c-atoms is fixed to be 
NLai = 5 + 5 = 200 and iV,=,^^,, = - = 60. 
Without impurities, as shown in Fig. 1(a), we have an LO 
state for the c-atoms with a period of A^: L — 16 along 
X and y-direction, which is close to the estimation value 
27r/g « 17.8 where q is the difference between the Fermi 
surface radius of spin-up and spin-down c-atoms. 

Interaction-induced localization-Despite the LO state 
in Fig. 1(a) breaks the translational symmetry of model 
(1), it preserves the C4 symmetry of the optical lattice 
and has a new translational symmetry which can be char- 
acterized by where q^^y — 27r/16. The density of the 
c-atoms has a two dimensional periodic distribution 
with local maximum (minimum) values corresponding to 
the peaks (saddle points) of A^. For weak Ucd, the wave 
functions for the impurities are extended and the impu- 
rity density nf is periodically distributed, preserving the 
C4 and ""^-translational" symmetry. As an example, we 
show nf for five impurity atoms in Fig. 1(b). However, 
once Ucd exceeds a critical value 14, the C4 and "~^- 
translational" symmetry is broken and nf is limited in 
several small real-space ranges each of which has a radius 
r(l3| centered around the saddle points of the LO super- 
fluid, as shown in Fig. 1(c), suggesting that the impurity 
atoms are localized. We find Vc is weakly dependent of 
Nd and the value of Vc can be found from the plot of 



the inverse participation number = Yii{nfY / (S^nf )^ 
which characterizes the degree of the localization|l3|. 
As illustrated in Fig. 2(a), for a single impurity atom, 
when Ucd is smaller than Vc, is almost zero while 
it gets significantly enhanced once Ucd > Vc- Accord- 
ingly, the radius r, labeled as blue circles in Fig. 2(a), 
becomes smaller with the increasing of P~^, characteriz- 
ing a stronger localization. Notice the step-like data lines 
of r in Fig. 2 is due to the lattice effect. From the inset 
of Fig. 2(a), we find that Vc increases monotonically with 
td, indicating that a stronger interaction is needed to get 
a lighter impurity atom localized. 

To explore the origin of the localization, we plot all the 
three parts of the ground state energy as the function of 
Ucd in Fig.2(b), in which Ec is the energy of c-atoms, 
Ed is the kinetic energy of the single impurity atom and 
Ecd — 'Yliii^cdnlnf) is the interaction energy between 
them. When Ucd < Vc, the energy of the c-atoms and 
the kinetic energy of the impurity atoms are almost un- 
changed while Ecd increases linearly with the increasing 
of Ucd- Once Ucd > Vc when the localization happens, 
Ec and Ed increases while Ecd increases in a smaller rate 
then decreases with Ucd- This indicates that the extended 
state is favored by Ec and Ed while the localized state is 
favored by the interaction energy and the localization is 
a direct consequence of the competition between them. 

Here the impurity atoms can get localized at the posi- 
tions where the potentials are distorted by themselves 
through the interaction Ucd- This can be interpreted 
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Figure 2. (Color online) (a): Plot of and the radius r as 
the function of Ud- Here td = 0.3. Inset: Plot of Vc as the 
function of td for model (1). (b): Plot of Ec, Ed and Ecd as 
the function of Ucd- (c) and (d): Plot of r and total energy 
Ed as the function of (c): j3 with td = 0.2, a = —2.8 and (d): 
td with Q = —3, P = 10.8 for model (4). In (d), tlo indicates 
the radius from model (1) with Ucd = 2.6. In (c) and (d), the 
lattice size for model (4) is 50 x 50 and the fitting parameter 
are A' = 10.6, B' — 2.0. Notice in all the figures the data 
lines correspond to the coordinates with the same color. 
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as a process of self- localization 3" 1^, which is in anal- 
ogy with the formation of small polarons in condensed 
matter physics where the electrons can get self-localized 
by strong electron-phonon interaction[l8|. For an intu- 
itive understanding of the self-localization of an impurity 
atom, we propose an effective model by assuming that the 
density of the c-atoms can be expanded as a function of 
nf, i.e. Ucdn"^ — Cn — anf + /3(f^f )^ + ... where C„ is a 
function of the density distribution of the c-atom in the 
absence of the impurity and does not change the nature 
of the self-localization qualitatively. By ignoring the ir- 
relevant terms and up to the second order of nf, we have 
the following effective model for the impurity atom: 



H, 



-td 



Pinff]dld,. (4) 



<h3> 



where a, /3 > 0. The most important feature of model 
(4) is that the interaction part between the impurity and 
LO state in model (1) is reduced to an self-consistently 
determined external potential. This strategy is in anal- 
ogy with the Kohn-Sham density functional theory |I19|, 
although different physical origin is presented there. The 
main physics of the self-localization can be grasped by 
this effective model by reducing several competing ener- 
gies to two parts: the first kinetic energy term and the 
second a- or "external potential" term, where the former 
favors an extended state and the latter favors a local- 
ized one. There exists a critical value ac at which the 
self-localization of the impurity atom occurs. A positive 
P term stabilizes the localized state with finite localiza- 
tion length A. To make an estimation on A, we assume 
that the total energy of one impurity atom from model 
(4) should be = At^/A^ - Ba/>^ + Cl3/X^ where 
A, B and C are integral constant [20|. By minimizing 
the total energy we have the optimal localization length 
A = ^J2C /B^J PI [a — A/Btd) which is proportional to 
21| . From Fig.2(c)-(d) we can see that the numerical 



results for r as the function of f3 and td are well fitted by 
our analytical estimation. With the increasing of r, the 
localization becomes weaker and the energy of the local- 
ized atom Ed is closer to the extended case with energy 
— 4trf. We also plot the variance of the radius from model 
(1) (labeled as r^o) as the function of in Fig. 2(d) and 
find its behavior is qualitatively consistent with the re- 
sults from our effective model. 

Paired state-When the impurity number is slightly 
larger than the saddle points of the LO state, they could 
form pairs. In Fig. 3(a), we plot the density distribution 
of nine impurity atoms. We can see that seven atoms 
are localized at different saddle points separately, while 
two atoms form a pair at one saddle point, which can 
be identified by their energy levels and wave functions. 
From Fig. 3(b) we can see that the paired states occupy 
the highest and lowest energy level, and their probability 
densities |^'i^9(i)p are also different from the unpaired 
states, as shown in the inset-plots of Fig. 3(b). Further 



calculations (not shown here) indicate that the paired 
state also forms on a lattice with bigger size, excluding 
the possibility of the finite size effect. 

The pair formed by two impurity atoms is also the 
consequence of the interactive process between the im- 
purities and LO superfluid. To see this, we make the 
ID-cut-plot of n° and nf for the paired state along its 
long axis. From Fig. 3(c) we can see that, different from 
two fermions trapped in a "hard" potential well, two im- 
purity atoms here modify the density distribution of c- 
atoms and create a double-well structure (labeled as red 
arrows) through the interaction Ucd- The inter- well tun- 
neling splits the energy into two levels which has the wave 
function of bonding and anti-bonding state, respectively. 
Due to their fermionic nature, two impurity atoms oc- 
cupy the bonding and anti-bonding state simultaneously 
and their wave functions ^Pi^g can be well fitted with 
{(f>L ± (/)i?)/-\/2 where 4'l,r is the wave function of the in- 
dividual well without tunneling, as exhibited in the inset 
of Fig. 3(c). Compared to the unpaired ones, the paired 
impurity atoms induce stronger local distortion to the 
LO superfluid, as shown in Fig. 3(d), which is unfavored 
by the c-atoms. Therefore the impurity atoms do not al- 
ways form paired states since they may not be supported 
by the total energy. 

Impact of mobile impurity atoms on the LO superfluid- 
As demonstrated in Fig.4(a)-(c) where Nd — 40,51,64, 
the LO states are significantly distorted. Although it is 
not unusual to find the symmetry reduction of the LO 
state caused by the mobile impurities here, a surpris- 
ing fact is that in all the three cases LO state maintain 
its original period to a certain degree or has a "quasi- 
"^-translational" symmetry. This phenomena identifies 
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Figure 3. (Color online) (a): 3D plot of nf for nine impurity 
atoms with td ~ 0.2 and Ucd ~ 5.0. (b): The energy levels 
of the impurity atoms in (a). Inset: |*I'i,5,9(j)p. (c): ID cut- 
plot of nf and nf for the paired states along the long axis. 
Inset: ^i(j) and ^'9(4). (d): Ai for the case in (a). Notice 
the positions of the impurity atoms are labeled as " x " . 
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an important difference between the mobile impurities 
and general fixed-impurities, which are randomly dis- 
tributed and always destroy the ""^-translational" sym- 
metry of the LO stateji^l- From Fig. 4(b) we can see 
that the symmetry of the LO superfluid with Nd = 51 
is lower than the case of Nd = 40. This is because 
the number of the impurity atoms Nd — 51 does not 
support a commensurate structure of impurity atoms 
with the LO state, resulting in more distortions for 
the LO superfluid. Compared to the former two cases, 
the pairing order parameters of the superfluid with 64 
impurity atoms show a quasi-one-dimensional-stripe-like 
pattern with non-crossing nodal lines. This indicates 
that the mobile impurity atoms can also affect the pat- 
tern of the LO superfluid besides the Zeeman field and 
temperature[2^ . 

Since both the polarized c-atoms and the impurity 
atoms are favored to stay at the saddle points, frustra- 
tions caused by their competition to the structure of the 
magnetization are inevitable, and thus novel distortions 
of the LO state can be expected for the cases of high im- 
purity densities. On the other hand, LO states with dif- 
ferent patterns caused by the impurity atoms (e.g., the 2D 
checkerboard and quasi- ID pattern) can be distinguished 
by detecting the spatial distribution of the population- 
imbalance experimentally. As shown in Fig.4(d)-(e), for 
a 2D-like LO state, although suppressed by the impu- 
rity atoms, the local maximums of the magnetization 
still have crossing points which correspond to the saddle 
points of Ai. For a quasi-lD LO state, the magnetiza- 
tion at the previous crossing points are totally destroyed 
and the polarized c-atoms are squeezed into non-crossing 
nodal lines, as shown in Fig. 4(f). 

When Ucd is increased, the impurity atoms tend to 
localize near the saddle points or nodal lines, and re- 
pel out the local c-atoms of the background to minimize 
the interaction energy. This is similar to the picture of 
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Figure 4. (Color online) (a)-(c): Contour-plot of Ai of the LO 
state with (a): Na = 40, Ucd = 2.5, (b) Na = 51, Ucd = 3.5 
and (c) Nd = 64, Ucd ~ 5.0. (d)-(f): Contour-plot of rrii with 
the same parameters with (a), (b) and (c), respectively. Here 
td = 0.02 and Uc = 3.7. 



phase separation. A pure phase separation where im- 
purities and superfluid occupy spatially different regions 
may be possible in cases with large impurity number and 
very large Ucd, which is beyond the scope of the present 
work. With and without impurities, the superfluid state 
is found to be LO state, whereas other states, like FF 
state and Sarma state, are never found stable in the pa- 
rameter space we studied. Although there is no long 
range order in a pure 2D infinite system at finite tem- 
peratures, our results should be observable at finite tem- 
peratures because the cold atoms are confined in a 2D 
optical lattice and thus may form a quasi- 2D system. In 
such a system the fluctuation effect is much suppressed 
and our mean fleld theory is expected to be valid. 

Proposal for experimental setep-Experimentally, 
Fermi-Fermi mixture was realized by using ^Li and 
'*''K[24] and the former was used to realize the possible 
FFLO state 0-il- 

Therefore two component ^Li atoms 
with spin-imbalance can be used as the background, 
i.e., the c-atoms, while one component ^''K atoms can 
be treated as the impurities. Considering is much 
heavier than ^Li, the critical value of Ucd for localization 
should be relatively small and our mean-field treatment 
is applicable. 
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Note added-Aitei completion of this work, we noticed 
a similar work in which the velocityof a mobile impurity 
was investigated in one dimension|251]. 
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